Impurity-induced states on the surface of 3D topological insulators 
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We calculate the modification of the local electronic structure caused by a local impurity on 
the surface of a 3D Topological Insulator. We find that the LDOS around the Dirac point of the 
electronic spectrum at the surface is significantly disrupted near the impurity by the creation of 
low-energy resonance state(s) - however, this is not sufficient to (locally) destroy the Dirac point. 
We also calculate the non-trivial spin textures created near the magnetic impurities and discover 
anisotropic RKKY coupling between them. 
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The Dirac spectrum of chiral excitations are real- 
ized in a wide range of materials including d-wave 
superconductors [JJ, grapheneJS], and semiconductors [3]. 
The Dirac spectrum brings in substantial similarities in 
electronic properties - like response to defects as well as 
low energy and low temperature properties. It is thus 
natural to combine these materials into a category of 
'Dirac materials'. A recent exciting realization of the 
Dirac spectrum is on the surface of 3D Strong Topological 
Insulators (STI)0H3- These materials have an ungapped 
spectrum at the surface while being fully gapped in the 
bulk. In addition, STIs are unique because the topol- 
ogy of their bulk band structure constrains their surface 
states to possess an odd number of Dirac nodes [7] ■ Sup- 
pressed backscattering inside the odd Dirac cone guaran- 
tees that the Dirac dispersion remain essentially unper- 
turbed for any perturbation to the Hamiltonian that pre- 
serves time reversal symmetry. This is a manifestation of 
the topological protection enjoyed by this kind of surface 
band crossing and makes these materials an attractive 
candidate for spintronics applications [5] as well as a pos- 
sible platform for topological quantum computation]!)]. 
In this work we focus on the stability of the STI surface 
nodes as well as the modification of surface states around 
local impurities. We consider local potential as well as 
magnetic impurities and calculate the modifications to 
the Local Density of States (LDOS) and the spin density 
near the impurity site. These should be accessible by 
STM measurements. We find the following, 
(i) There is substantial modification of the LDOS near 
the impurity site for both the nonmagnetic (time rever- 
sal preserving) and magnetic impurities (time reversal 
breaking), especially when impurity scattering is strong 
(unitary). Near the potential/magnetic impurity, a sin- 
gle/a pair of low energy resonances form near the Dirac 
point (Figs. [lj[2|- These become very sharp and their 
energies f2 — > as the impurity strength ((3j) |Z7| — oo : 



The scalar impurity resonance is doubly degenerate due 
to Kramers' theorem [Tl [TUI ITT] . The magnetic impurity 
breaks time reversal symmetry and splits the low energy 
impurity resonance into two spin-polarized resonances on 
either side of the Dirac point (Fig. [I]), 
(ii) Modification of the LDOS vanishes quickly for en- 
ergies less in magnitude than the resonance energy (ap- 
proaching the Dirac point) for both magnetic and non- 
magnetic impurities. Thus, modifications to the low en- 
ergy LDOS does not provide us with a signature of any 
incipient gap in the spectrum for both potential and mag- 
netic impurities. For r >> l/ui, these decay as 1/r 2 . 
(in) In addition to LDOS modifications, magnetic scat- 
tering produces non-trivial spin textures near the impu- 
rity site (Figure [3]) that can be imaged with a magnetic 
force microscope or spin-resolved STM. These non-trivial 
spin textures lead to the propagation of unconventional 
antiferromagnetic (AF) RKKY coupling between mag- 
netic impurities, when they are polarized along the line 
joining them and when the chemical potential is close to 
the Dirac point. When the spins are perpendicular to 
the line joining them, they interact strongly and ferro- 
magnetically (FM). The Dzyaloshinskii-Moriya (DM) in- 
teraction between the spins vanishes at the Dirac point. 
We thus conclude that random magnetic impurities will 
tend to align parallel to the normal to the STI surface. 

We will model the STI surface states as a single species 
of non-interacting 2-D Dirac quasiparticles[12] with a 
high energy band cutoff W. We shall work in units of 
W, h and Vp (the Fermi velocity). The Hamiltonian be- 
comes: 

U = ct-p (2) 

where er/2 is the actual spin of the electron (or related 
by a rotation about z). We shall consider local impurities 
of the potential and classical types respectively: 

V pat = U 15(f), K„ ag = US ■ <r6(r) (3) 



(1) 



For the magnetic case we have assumed a local Heisen- 
berg exchange J between the band electrons and the im- 
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FIG. 1: (Color Online) LDOS plots showing the low energy resonance(s) near (A) the scalar/potential impurity, (B) the z- 
polarized and (C,D) the x-polarized magnetic impurities. (C) and (D) show the a;-spin projected LDOS, at a point on the x 
and j/-axis respectively. Note from (C) that on the x-axis, the negative energy states have excess states with spins parallel to 
the x-polarized impurity. In all these cases, U = 100, r = 20. In the system of units used above, h, vp and W are unity. 



purity spin S, whose direction is given by the unit vector 
S. Thus, U = JS/2 in F mag . 

To address the effect of impurity scattering we use the 
T-matrix technique pQ. The T-matrix is defined via: 



T(w) = V + VG' et (u)T(uj) 



(4) 



where Gg et is the retarded Green's function for the 
impurity-free material and w is the energy. For uj <C 1 
and p ^> 1/W (p = r — r'), it has the following form: 

(r |GTMI O = ^ [fo(", P)l + /iK p){a ■ p)] (5) 
where 

/o(w,p) = s( w )r - U O 0,fi(u,p) = iYx+s{oj)J x 9 (6) 

and \ui\p is the argument of the Bessel functions J /i and 
y 0/1 . Also, s(-) ee sgn(-) and6> ee 8(1-|w|). We shall also 
require the unperturbed on-site Green's function valid for 
short distances < 1: 

GqH = (0|G^H|0) =-(g {u))+i 9l (u))I, where 



9o M = 



- 1 



^e(i-M) (?) 



In pi, we have used a local form for the impurity po- 
tential (t-| t^jr*') = V S(r)5(r'), where V is a 2 x 2 matrix 
in spin-space. The T-matrix also becomes (r\T\r') = 
T S(r)S(r'), with T satisfying the following equation 



T = V + VG T = (I - VG )~ 



(8) 



From the algebraic relations involving (J3j) , Q and p]), 
we analytically calculate the T-matrix, the full Green's 
function G rot , 

G-» = GTM + Gr(o;)f (9) 
the full (spin-unresolved) LDOS, 



p(r,w) 



-— ImTr < r 

7T 



G Iot (u;) 



(10) 



the local density of spin up/down states (in direction //), 



p± (r, w) = ImTr ( r 



and the energy-resolved spin density averages: 



s(r,uj) 



--ImTHr 

7T 



<S~(«)f 



(11) 



(12) 



For the scalar and magnetic impurity cases, we find that 
the additional GTG = SG pieces in the Green's function 
([9]) evaluate respectively to (using g = g\ + ig%) 

5G pot = ^ft_^ (13) 



and 
SG 



Ulu 2 



16 l + Ug 



■[-2if f l( T-(Sxry 



mag 16(1 - U 2 g 2 ) 
(f 2 + ! 2 )a ■ S - 2/ x V • r)(S ■ f) - Ug(f 2 - f 2 )] (14) 
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FIG. 2: (Color Online) Low energy LDOS near the (A) scalar 
impurity and (B) the a;-polarized magnetic impurity for U = 
80. h, v and W have been set to unity. 



As shown in Figures ([!]) and ([2|, for both the mag- 
netic and non-magnetic cases we obtain low energy res- 
onance^) in the LDOS (arising from the minima of the 
denominators in (13 1 and ( 14 ) ) that approach the Dirac 
point for large impurity strengths according to ([I]) [TlllOj. 
These resonances become sharper as they approach the 
Dirac point with increasing potential strength and while 
doing so, also increase in amplitude relative to the un- 
perturbed LDOS. For r ^> and u<Cl, the strength 
of LDOS modulations diminish with distance as l/r 2 |13|. 

Topological stability of the surface Dirac spectrum in 
TIs is often discussed as a crucial property of these ma- 
terials. An important question in this context is whether 
the appearance of these low-energy resonances is related 
to the local creation of a gap at /destruction of the Dirac 
point. Naive scaling analysis tells us that the potential 
strength U has a dimension of —1 (same as length) near 
the fixed point corresponding to As we approach the 
Dirac point, we should thus see the effects of the impurity 
become negligible. Indeed, we find that if we move from 
the resonances to the Dirac point, the density of states 
gradually settles down to the impurity-free value. We 
cannot, therefore, find signatures of gap-opening at the 
Dirac point at the stage of one-impurity scattering. We 
also note here that the appearance of these resonances 
at the Dirac point is a consequence of the band cutoff 
being symmetric on the particle and hole sides - in re- 
alistic materials [1] the band structure is asymmetric and 



depending on the degree of asymmetry, these resonances 
may appear at other region(s) of the bandspQ. 

In addition to the impurity resonances at small energy 
we find new states that lie outside the effective band 
edges - a consequence of using a hard cutoff. These 
true bound/anti-bound states are located at the posi- 
tive/negative side for positive/negative sign of a scalar 
impurity potential U. For a magnetic impurity they are 
located on both sides outside the effective band edges. 
For large \U\, these are located approximately at a dis- 
tance U from the Dirac point, while as \U\ — > they 
approach the band edge as e \ u \ . In real STI SSs, these 
may well be located at the same energy as the bulk bands, 
will hybridize with them and delocalize into the bulk. 

Near a magnetic impurity, entanglement of the electron 
spin and momentum lead to the creation of spin textures, 
as shown in Figure [3] The energy-resolved spin average 
is found to be: 



s{r,Lu) (15) 

' 2i fofiS xr - (/ 2 + fl)S + 2f 2 r(S ■ r ) s 
16tt(1 - U 2 g 2 ) 



= Uuj z Im 



The first term in ( 15 ) gives rise to a DM interaction be- 



tween two impurity spins. When the chemical potential 
/j, is at the Dirac point, considering only the perturbative 
result (obtained cheaply by putting g — ¥ in the above 
expression) , the strength of this interaction becomes zero 



Re(/ /i)w 2 dw 
1 f°° d 



r3 Jq dx 



(J y )O(r - x)x z dx « 



(16) 



For a finite chemical potential -C 1, the amplitude of 
the DM interaction becomes 

r / " R e (/ o/l )c 2 rfc = ^^ [ M Mfoh)<Sd»> 



8tt 



Stt 



(17) 



At large distances, the amplitude of this interaction de- 
cays as ~ Ufi/r 2 . 



The second term in ( 15 ) leads to FM RKKY interac- 



tions when fi = 0, in the perturbative approximation. 
The corresponding spin component is 



US 

167T 



Im(/ 2 + /jVdw 



U 



327T7- 3 



(18) 



Finally, the third term in (15) leads to AF RKKY 



interaction between impurity spin components pointing 
along the line joining the impurities, when fi = 0. In 
the perturbative limit, the corresponding induced spin 
component is: 

Uf{S - f) f° Im(/> 2 ^ = ^r(S.r) (19) 



NTT 



4 





(A) 






(B) 




























• . / / / \ \ \ - • - 












: : : ;l 




Y , 




• ■ ■ » v\ f , , . 


• ' J 






o T 




. . ■ . \ \ 




{!:::: 




.,,//< \ \ > . ■ 


. . ■ . \ \> 








..,11 \^/ ( 1 , . . 










- > \ \ \ •v — *- / 1 I ' • ■ 








-20 





FIG. 3: (Color Online) Spin textures near spin impurities 
(U = 20, E = ft = -0.08) when the impurity is (A) z- 
polarized and when it is (B) x-polarized. The component 
in the xy plane is denoted by a vector while the background 
shade gives the sign of s z (r, E) (clear = positive) . The arrows 
are normalized to the longest field-of-view total spin length 
in (A) and xy spin length in (B), indicating respectively the 
sign of the z polarization and the anisotropic x polarization 
around the impurities in accordance with ( |15[ ) (these mediate 
anisotropic RKKY interactions between two impurity spins). 



From these two expressions we conclude that at ji — 
the interaction energy gained when two impurity spins 
point at/away from each other is half that gained when 
they are aligned parallel to each other and perpendicular 
to the line joining them. Thus, when many impurities 
are present (and /z = 0), they will tend to point in the 
common direction where all gain the FM interaction en- 
ergy - along the z-direction, normal to the surface. This 
kind of FM ordering will be conducive to opening a gap 
in the STI surface state spectrum. 

We would like to note here that the foregoing results 
are not obvious when observing the energy-resolved spin 
densities at low energies to — > 0— , because of the low den- 
sity of states there. Naively, one would have expected 
the low energy long wavelength features to determine 
the r — > oo spin textures, but the low energy spin tex- 
tures predict, incorrectly, antiferromagnetic RKKY in- 
teractions with short distance ferromagnetic contribu- 
tions arising from non-perturbative effects. We would 
also like to note here that we assumed a smooth cutoff 
when adding up the spin textures at different energies 
to eliminate cutoff-dependence [13]. We have used mul- 
tiplicative functions like e - ' 7 '"'' and (1 — e~ vu )/(rjuj 2 ) 
(having different characters as \lu\ — > oo) in the energy 

integrals and then taken the limit r\ — > OH both these 

procedures gave the same limit [16 . 

When we consider the full nonperturbative spin aver- 



ergy resonance(s). However, in the asymptotic approach 
to the Dirac point, the linear DOS is preserved, consis- 
tent with the negative scaling dimension of the impurity 
strength. Thus, the gap-opening mechanism for mag- 
netic impurities is not evident at this stage of analysis. 
We also find that the induction of non-trivial spin tex- 
tures near magnetic impurities leads to the mediation of 
antiferromagnetic RKKY coupling between impurity spin 
components parallel to the lines joining them, especially 
if the chemical potential is at the Dirac point (in which 
case the interaction does not oscillate in sign). The spin 
components perpendicular to the line joining the impuri- 
ties, however, exhibit strong FM interaction. While there 
is, in general, a DM component in the spin interactions, 
it vanishes at the Dirac point. Finally, when many mag- 
netic impurities are present (and fi = 0), they will tend 
to align ferromagnetically normal to the surface. 

We are grateful to D. Abanin, D. Basov, Z. Hasan, 
H. Manoharan, N. Nagaosa, T. Wehling and Y. Xia for 
useful discussions - especially to D. Abanin for drawing 
our attention to the method used in [TS] to calculate the 
RKKY interactions. This work was supported by the US 
DOE thorough BES and LDRD funding and by Univer- 
sity of California UCOP program T027-09. 



age obtained by integrating ( 15 ) numerically, the afore- 
mentioned behaviors seems to hold qualitatively if we 
look beyond the 'ringing' introduced by a sharp cutoff. 

In summary, we find that local impurities can strongly 
disrupt the structure near the Dirac node of 2-D surface 
states in 3-D topological insulators by forming low en- 
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